Electronic pressure on ferromagnetic domain wall 
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' The scattering of the eletron by a domain wall in a nano-wire is studied perturbatively to the 

04 , lowest order. The correction to the thermodaynamic potential of the electron system due to the 

^ ' scattering is calculated from the phase shift. The wall profile is determined by taking account of 

Q I this correction, and the result indicates that the wall in a ferromagnet with small exchange coupling 

. can be squeezed to be very thin to lower the electron energy. 

' ' Recently the magnetoresistance in various magnetic nano-structures such as multilayers, wires and contacts has 
, been attracting growing interest. The magnetoresistance effects in such systems are induced by the spatial variation 
' ' of the localized magnetic moments. Hence, to grasp the physics of the resistance due to a domain wall scattering is 
fundamental to an understanding of the mechanism of the magnetoresistance. The theoretical study of this problem 
' goes back to more than 25 years ago by Cabrera and Falicov [Q. The scattering mechanism by the wall is due 
' to the non-adiabatic deviation of the electron spin from the local magnetization. In the case of a thick wall with 
. thickness A ^ k'^^ {kp being the Fermi wave length), the electron spin can adiabatically follow the spatial rotation 
t/o ' of the magnetization as it traverses the wall, resulting in a negligibly small reflection. Cabrera and Falicov thus 
1^ concluded that the resistivity contribution due to domain walls is important only in a limit of fc^A — > 0. This limit 
H ■ was unrealistic at that time, where 3d metals in the bulk were considered mostly. However, recent technology has 
I t ' put the problem in the limelight, firstly by the direct measurements of domain wall resistance [^-Q and secondly by 
fabrication of a very small domain wall in nano-scale contacts |^,^. Motivated by the recent experiments detailed 
^ ' theoretical studies including the quantum correction in the dirty case 0,^, spin-dependent lifetime etc. have 
been made [| 10 -1^]. These investigations revealed that the problem was not as simple as had been believed. In fact 
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theories predict that even the sign of the domain wall contribution to the resistivity depends on temperature, system 
size, amount of disorder as well as the material, and experimental results indeed show such diversity. A clarifying 
^ ■ explanation of the whole behavior observed in wires has not given so far. In nano-scale contacts, the situation seems 
to be simpler. There the wall is believed to be constrained in a small contact region of an atomic scale due to a rapid 
change of the contact shape, and thus adiabaticity does not hold, resulting in strong scatterings |l^,^,|l^. A recent 
experiment on Ni contacts revealed a magnetoresistance of about 300% at room temperature ^ , which is significantly 
. large compared with other giant magnetoresistive materials such as multilayers. The result has been interpreted as 
f — ■ due to a strong electron scattering by a domain wall, which is trapped in the narrow contact region in the absence of 
T-H . the magnetic field, and which is expelled by applying a magnetic field 

' The strong scattering by the wall in nano magnets suggests that the electronic state is modified. This possibility is 
studied by looking at the phase shift, = Arg i, t being the transmission coefficient due to the domain wall scattering. 
This phase shift is related to the total energy of the electron system through the equality D = 111,0 , where D 
is the density of states and e is the energy. Thus the stability of the electron system in the presence of the domain 
. wall can be discussed in term of the phase shift. 

' The relation between the scattering and the stability of the electronic state has been studied in the non-magnetic 
I , case. Experimentally it was shown in a Au contact that the tensile force of the contact as the contact is stretched 
'■^ ' is related with the conductance and becomes maximum when the conductance shows a plateau with a quantized 
^ ', value of integer times 2e^//i This result indicates that the electron system is stabilized when the conductance 
O I takes a quantized value. Such behavior has been explained considering the scattering by the potential energy in the 
constriction during the prolongation and taking account of the energy shift of the electron system, which is calculated 
from the scattering phase shift by use of the above relation ||l^ . 

In the case of ferromagnetic wires and contacts, a domain wall may affect the stability of the electrons and thus the 
wall structure can be modified to lower the electron energy. The essential parameter of the wall which governs the 
' transmission coefficient is the thickness, A. In this paper the thermodynamic potential of the electrons is calculated 
as a function of A. This thickness is determined to minimize the total energy of the wall, i.e., the magnetostatic 
energy plus the electron energy. An interesting point is that in 3d metals the energy scale of the electron, ef, is much 
larger than that of the magnetization, the exchange (J) and anisotropy (K) energies. Thus even a small effect of the 
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scattering of the electron might results in an energy shift comparable to the magnetostatic energy. In fact it will turn 
out that this is the case in a ferromagnet with small J, and the wall would be compressed to be very thin, A ~ kp^. 

Considering a narrow wire, the system is treated as in one dimension. The calculation is at zero temperature. 
The electron channels considered here are those which contribute to the transport and the quantized behavior of the 
conductance is not discussed here. We consider the exchange interaction between the electron spin and the local spin 
and the Hamiltonian is 



kcr 



g / dxS{x){c^ (Tc), 



(1) 



where = fi^k^ /2m — ep {ep being the Fermi energy). The spin index is denoted by cr = ± and cr is the Pauli 
matrix (spin indices are suppressed in the second term). Here, the local spin S has a spatial dependence of a domain 
wall. In terms of a polar coordinate [Sz = ScosO) the wall is represented as cosO 1 for x ^ —A and cos6' — 1 
for a; ^ A, A being the thickness of the wall. To proceed we carry out a local gauge transformation in the electron 
spin space so that the z-axis is chosen to be along the local direction of spin ~§ ||igt] . The transformation is written as 
Oct = crcos (I) Co- — i sin (I) c_o-, where the electron in the new frame is denoted by a^. The Hamiltonian is modified 
to be if = Hq + i/int where Hq = ^ ^kaa^fjO-ka and €ka = Cfe — crAis the energy of uniformly polarized electron with 



the exchange splitting A 
given by 



g\S\. The term iJjnt describes the interaction between the electron and the wall and is 
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where 



Aq = I dxe-*«^V6l, 



(3) 



is a domain- wall form factor and L is the system size. Note that V0 — > for |a;| — > oo and so the interaction is localized 
near the wall. The spin-conserving process in the rotated frame as the electron goes through the wall corresponds to 
the adiabatic process, where the electron spin follows the rotation of the magnetization, and the spin flip corresponds 
to the deviation from adiabaticity. 

The interaction Hi^t leads to a scattering of the electron and thus contributes to a resistivity due to the wall as has 
been discussed before based on a linear response theory 0,^. Here we study the phase shift due to the scattering, 
calculating the T-matrix perturbatively with respect to Hi^t ■ The wave function corresponding to the incoming state 
with the energy E and the spin a, which we denote as {a' is the final state spin), is calculated as 



^ E 



'IK X rp 
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where G%^, {x — x') 



is the free Green function and Tk', 



E< k'a'\T\k^a >, T = ^int i_gu 



The transmission and reflection coefficients are defined from the asymptotic behavior as 
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The result is 
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By use of these coefficients, the resistance due to the wall is calculated based on the Landauer's formula (with 
two-terminals) as = R — Rq, where 



R-' = {eyh)J2{\to 



(7) 
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and i?o = {h/2e'^) is the resistance without the wall. 

The diagonal (spin a) component of the transmission coefficient for the incoming wave vector of is calculated as 



t(T(7 {k(7 ) 



1 



2 3fc^ + fc^^ 
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(8) 



The diagonal component of the phase shift, which is defined as 0cr(e) = Imlnfo.c,(fc) (e = (?i^fc^/2m)), is obtained as 

,4fc2 + 2kq-2aA 



0<x(e) 
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2kq 
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where A = {kp^ — kp_)/2 and P denotes taking the principal value. In the adiabatic limit of large kpX, the phase 
shifts of two spin channels acquire Berry's geometrical phases |2^. The phase shift difference of the two channel 
become tt since the solid angle of the rotation of local magnetization is tt. 

The phase shift contributes to the change of the density of state of the electron according to the relation 
SDa-{e) — {l/iT){d(j)a{f-)/ de) Thus the thermodynamic potential of the electron system (defined as f^o = 

— ^ In Tre^'^^^"'^^^)) is shifted due to the scattering by an amount of 
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where Hq = ^ ^Fa4'a{^ = 0) is an irrelevant constant (independent of domain wall parameter), which we neglect. 
Thus we obtain 
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where fka = 0{ep — eka) is the Fermi distribution function at zero temperature {0(x) is the step function). We note 
that this expression of SQc is also obtained directly by evaluating the diagrams described in Fig. || j2^. In fact these 
contributions are calculated as 
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which is equal to (O). 

The expression ^1^) reduces after fc-summation to 
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To proceed, we approximate the profile of the wall to simplify the calculation. In the case of an ideal wall determined 
solely by the exchange and uniaxial anisotropy, the profile is cosfl = tanh(a:/A) and so Aq = tt/ cosh(7rgA/2). By 
noting that the essential feature of the wall is that the form factor Aq is large only for \q\ < A~^, we here replace 
Aq = 2 for \q\ < A^^ and Aq — Q for \q\ > A^^. (We have checked that this approximation does not change the result 
very much.) By this approximation the energy shift is calculated as 



6n^ 



where 



-epefikpX), 
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f{x)^ — LT + -(l-4a;^)ln 



2a; + 1 



2x- 1 



(14) 



(15) 



and 
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The minus sign of 5^c indicates that the electron system is stabihzed in the case of a thiner wall {Xkp — 1) which 
causes a strong scattering. This is because the two spin channels are mixed by the domain wall. The function 
f{x) grows at smaller x as plotted in Fig. ^j. The result ( p^ ) of the perturbative calculation is not justified at 
fci;'A < 1 . We thus estimated the function f{x) non-perturbatively by calculating the transmission amplitude by the 
recursive Green's function method . We have used discretized form of the Hamiltonian Eq. (|l|) with a tight-binding 
approximation. The result of function / is potted in Fig. || as dashed lines. The function / turns out to depend on 
C only very little in the region physically meaningful (i.e., kpX ^ 1), which is consistent with the perturbative result 

(p^. It is seen that the perturbative result (solid line) gives a good approximation to the full result. 

In the case of an insulating magnet the structure of the domain wall is determined solely by minimizing the 
magnetostatic energy. The magnetostatic energy is written as 

i^M=/^r^^(V.)^ + ^sin^.V (17) 



a V 2 2 

where J and K are the exchange and anisotropy energies per site, respectively and a is the lattice spacing (we assume 
that the anisotropy energy K includes the contribution from the dipole energy). In the ideal case the wall thickness 
is given by Aq ^ \J J / Ka. Here we consider the case where the wall is thinner compared with the bulk case, A <C Aq, 
where the effect of the electron can be important. We assume that the cross-section of the wire is only weakly 
varying in the length scale we consider, and thus neglect the effect from the wire geometry discussed in ref. The 



magnetostatic energy of the wall is thus dominated by the exchange term of eq. (17). By approximating the wall as 



V9 ~ 7r/(2A) for |a;| < A and V9 ~ otherwise, the magnetostatic energy is obtained as 

E^^ Jdx'I^m'^^^^J. (18) 

In the metallic case the total energy due to the wall is given by the sum of the electron part and the magnetic part 
{E^ = Sn, + Em): 

E^^—J-eFefikpX). (19) 

It is seen that the electron energy and magnetostatic energy compete with each other. The behavior of E^ is plotted 
in Fig. ^. Only the region of very small A is plotted since for larger A the functions vary monotonically. 
The parameter which determines the relative strength of the two energies is 

r^^O.lx^, (20) 

where a factor of 0.1 is to take into account the value of the function f{x) at a; ~ 1. The value of ij exceeding 1 
indicates that the wall energy can be lowered if the wall becomes very thin. In 3d transition metals such as Ni, 
ep ~ 4eV, J ~ O.leV and ( ^ 0.3, and thus 77 ~ 0.15 in the case of d-band. Therefore the electron part is not strong 
enough to modify the domain wall structure. As is seen from the figure, in the case of ^ = 0.3, if J is as small as 
^ O.OOlei?, the wall is compressed to become an atomic size, kpX ~ 1. In a strong ferromagnet with ( = 0.8 the 
crossover value below which the electron energy wins occurs around J ^ O.lGe^^. 

If J is smaller than the threshold value, the wall is very thin and so leads to an almost perfect reflection \t\ ^ for 
channels below the threshold and \t\ ~ 1 for those above. In other words, the wall is self-organized to give rise to a 
quantized value of the conductance, just as the wire shape does during the prolongation in non-magnetic wires p8[ . 
In recent experiment of Ni contacts clear quantization of the conductance in unit of 2e^//i has been observed under 
zero and small magnetic field [ p3| . A possible existence of a domain wall has been suggested there, but naively domain 
wall in nano-contact is expected to lead to a deviation from the conductance quantization since it gives rise generally 
to a transmission probability of T 7^ 1. The present mechanism might be relevant to the clear quantization observed, 
since suppression of the effective exchange coupling J^g is expected in nano-contacts due to low-dimensionality. (The 
suppression of Jes is indicated in experiments on films as a decrease of ferromagnetic transition temperature (Tc) as 
the thickness is reduced p^.) 

In ref. [^, the importance of the wire shape on the wall structure has been demonstrated. This effect is due to 
the shape anisotropy and so the relevant energy scale is J 01 K. Hence the compression of the wall by the electron 
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system would be more important in ferromagnetic wire with small J or large ^. We have calculated in one-dimension 
but the present result would qualitatively hold in the bulk case as well. 

Suitable system for experiments would be a low Tc ferro- or ferri- magnetic materials. The present study indicates 
that in the metallic case there can be a crossover value of J, below which A becomes of order of lattice spacing. From 
the point of view of controlling the electron density, experiments on magnetic semi-conductors might be interesting. 

To summarize, the energy shift of the electron system by the scattering by a domain wall has been calculated. The 
result indicates that the wall lowers the energy of the electron system and the effect is larger for a thiner wall than 
a thick one. This contribution thus competes with the magnetic exchange energy, which favors a thick wall. In 3d 
metals in the bulk the effect turns out not to be strong enough to affect the wall structure, but it can dominate in the 
case with smaller J or a larger splitting or in nano-wires, in which cases the wall is squeezed by the electron system 
to be about an atomic scale. 
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FIG. 1. The lowest order contribution of the wall to the electron energy. The solid and wavy lines denote the electron and 
the wall, respectively. 

FIG. 2. The function f{x). Solid line is the result of a perturbative calculation ((eq. (|l5|)) and dashed lines are the 
non-perturbative (numerical) results for ( = 0.3 and ( = 0.8. It is seen that / does not depend much on (, and the perturbative 
approximation is fairly good. The region kpX < 1, which is not physically meaningful, is plotted by dotted lines. 
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FIG. 3. The total energy due to a domain wall, E^/tp, as a function of thickness, A. Solid and dashed lines axe for = 0.3 
and 0.8, respectively, and three lines correspond to J/ep = 0.02, 0.01 and 0.001 from the top to the bottom. It is seen that 
the energy of the electron system is lowered for a very thin wall, /cfA ~ 1, if J /ep ^ 0.001 {J/^f ^ 0.02) for C, = 0.3 (C = 0.8), 
respectively. 
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